Large time off-equilibrium dynamics of a manifold in a random potential 
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We study the out of equilibrium dynamics of an elastic manifold in a random potential using 
mean-field theory. We find two asymptotic time regimes: (i) stationary dynamics, (ii) slow aging 
dynamics with violation of equilibrium theorems. We obtain an analytical solution valid for all large 
^ ■ times with universal scalings of two-time quantities with space. A non-analytic scaling function 

' crosses over to ultrametricity when the correlations become long-range. We propose procedures to 

I , test numerically or experimentally the extent to which this scenario holds for a given system. 

Oh' 

The dynamics of an elastic manifold in a quenched 
On| ' random potential is relevant for a large number of ex- 
perimental systems. Examples are flux lattices in high- 
Tc superconductors 0, interfaces in random fields 
charge density waves B , dislocation in disordered solids 
' pL surface growth on disordered substrates jffl. The 
. competition between elasticity and disorder produces a 
0^ ' 'glass' state that may exhibit pinning, slow dynamics and 
non linear macroscopic response (e.g, leading to zero lin- 
ear resistivity in superconductors Q). While there is 
a phenomenological picture based on scaling arguments 
(droplets) ||] no satisfactory analytical approach is avail- 
able at present for the low-temperature dynamics. 

The statics of a d-dimensional elastic manifold em- 
bedded in a N dimensional space in the presence of 
a random potential was studied by Mezard and Parisi 
£3 ■ (MP) who applied a replica variational Gaussian approx- 
imation (Hartrec) for N finite which becomes exact at 
• i-h , N — oo ||. The replica symmetry breaking (RSB) solu- 
' tion captures some of the essential physics in finite N di- 
mension, such as sample-to-sample susceptibility fluctu- 
ations and predicts non trivial (Flory-like) rough- 
ness exponent C- It allows for a theory of the statics 
of the vortex glass state in superconductors relevant to 
experiments such as flux decorations. Other analytical 
approaches are based on RG methods Jll|-p^| and it is as 
yet unclear whether they capture all the physics Q,^5[. 
Despite the obvious interest of the static approach, it 
applies by construction to equilibrium (Gibbs measure) 
properties, which may not hold for experimental times in 
a glassy system Q . 

In this Letter we study, also within the Hartree ap- 
proximation, the dynamics of this problem starting from 
a random configuration as in a temperature quench. We 
find that at low enough temperature there is an aging 
regime and that the system never reaches equilibrium. 
The correlation functions then depend not only on time 
differences but also on the waiting time after the quench. 
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We obtain the two-time scaling with explicit space depen- 
dence, a new feature with respect to the mean-field ana- 
lytical results for glassy systems obtained so far [jl?] 5oJ . 
Details will be presented elsewhere . 

This treatment is exact for TV = oo. At finite TV it holds 
a priori only within the Gaussian variational ansatz. 
Whether the properties hold qualitatively for true finite 
TV dimensional models cannot be answered analytically 
at present. The good qualitative agreement Q of the 
MF dynamical analytical solution with experiments on 
spin glasses suggests that our results for TV = oo may be 
relevant for some systems related to the present model. 

The main purpose of this Letter is to suggest, on the 
basis of the exact solution for infinite TV, definite pre- 
dictions for the off-equilibrium dynamics which can be 
checked in numerical simulations and experiments. Our 
results provide a basis for a finite-TV analysis. Remark- 
ably, despite the system being out of equilibrium, some 
of the results of the replica calculation by MP,e.g £, are 
shown to carry through to the dynamics, albeit with a dif- 
ferent interpretation in terms of directly observable time- 
dependent physical quantities. 

The model of a manifold of internal dimension d em- 
bedded in a random medium of dimension TV is de- 
scribed, in terms of an TV component displacement field 
<fi, <fi — (4>ii $2, ■ ■ ■ , <Pn)> by the Hamiltonian H] 



H = d d x 



l -(y<t>{ x )f + v{<t>{x), x ) + ^<t> 2 

fi is a mass, which effectively constraints the manifold to 
fluctuate in a restricted volume of the embedding space. 
V is a Gaussian random potential with correlations 



V{t(>, x)V(cf>', x 1 ) = -N8 d (x - x') V ((</> - fi f/N) 
We consider the Langevin dynamics 
d(j>(x,t) _ 8H 



dt 5<f>(x,t) 



T]{x,t) 



with (t) Q (x, t)T) P (x', t')) = 2T S aP S d (x - x')S(t - t'). 
We let the system evolve from a spatially translationally- 
invariant (STI) configuration at t = 0. The system re- 
mains STI at subsequent times. We study the dynamics 
exactly for TV — > oo and, for TV finite, within the dynami- 
cal Hartree approximation which amounts to substituting 
V by an effective V fp| |2fj|. The 'equilibrium' dynamics 
(d la Sompolinsky [P2|) was studied by Kinzelbach and 
Horner for general d |23| . The case d = has been stud- 
ied in Ref. |Q and an analytical solution at large time 
was given in Ref. |2Q1 . 

The quantities of interest in the large-time off - equi- 
librium dynamics are the two-time correlation C xx > (t, t') 
— 1 /N((p(x, t) ■ 4>(x' , t') ) and the response R xx '(t,t') = 
l/Nd(<j>(x,t))/6f(x',t')\ f=0 where f(x',t') is a small 
perturbation applied at the space point x' at time t'. We 
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also define the mean squared displacement D xx i(t,t r ) = 
1 /N{((j> {x, t) ~ 4>(x',t')y) and the correlation B xx . (t,t') 
= l/N{{cj>{x,t)-<t){x,t')){ci>[x',t)-ci>{x',t'))). The 
brackets and overline represent the average over the ther- 
mal noise and the quenched disorder, respectively. The 
Fourier-transform Bk and R k w.r.t. the space difference 
x — x' are used below, and we denote with tilde equal 
space (x = 0) two-time functions, B = j k Bk- 

A common choice is V(z) = (9 + z) 1_7 /(2(l — 7)). 
'Short-range' (SR) correlations correspond to 7(1 — 
d/2) > 1 and 'long-range' (LR) correlations to 7(1 — 
d/2) < 1. The static solution II is characterized by two 
exponents: in the LR case, an ultrametric ansatz gives 
D x * = (4>(x) — (f)(0)) 2 ~ x 2 ^ with a roughness exponent 
Clr = (2 — d/2)/(l + 7) and a free-energy fluctuation 
exponent 9. In the SR case a one step RSB ansatz gives 
Csr = (2 — d)/2. Statistical rotational symmetry imposes 
in general 9 = 2( + d — 2. The d = 2 Sine Gordon model 
is marginal and solved by a one step RSB. 

Aging. Let us describe the picture that emerges for 
large times and fi > 0. Equal-time quantities reach their 
asymptotic values which do not necessarily coincide with 
the equilibrium ones ]T^ , |2l[ | . For two-time quantities we 
consider two different regimes of the times. 

(i) After a large waiting-time t w , Bk(r + t w ,t w ) first 
grows with r in a manner independent of t w , from up to 
the Edwards- Anderson parameter for the mode k defined 
as b\ = lim r ^ 00 lim^^oo Bk(r + t Wl t w ). In this time 
regime the displacement is time-translation invariant 
(TTI); we denote it 6f(r) = lim^-^ B k (r+t w , t w ). The 
response r^(r) = lfm^ _,oo Rk(r + t w ,t w ) satisfies fluc- 
tuation dissipation (FDT) rf (r) = 1/(2T) d T b£ (r) 9(t). 
The 'FDT-regime' is very much like an equilibration in 
a state, the manifold looks pinned with an effective mass 
M, since b\ = (2T)/(k 2 + M) (for fi -> 0). 

(ii) However, for all t w and sufficiently large r, 
Bk(r +t w ,t w ) continues to grow beyond b\ up to b% = 
limr—xx) Bk(r + t w ,t w ). The growth of B^ now depends 
on t w : the larger t w the slower the motion of the system, 
it ages. Thus in this model pinning is a gradual pro- 
cess, the older the system the more pinned it is but it is 
not pinned forever. The aging time-regime is thus defined 
(for t w — > 00) as the times r such that Bk(T+t w , t w ) > b\. 
We denote by &&(t + t w ,t w ) and rfe(r + t w ,t w ) the dis- 
placement and response in the aging regime, where both 
TTI and FDT are violated. 

As regards the measurements of noise and susceptibil- 
ity (i) corresponds to high frequencies while (ii) corre- 
sponds to low frequencies (scaling with the waiting time) 
such that noise and susceptibility depend on t w . In a 
domain-growth process (i) corresponds to the fast ther- 
mal fluctuations of the spins around their mean magneti- 
sation while (ii) corresponds to the actual growth of the 
domains. An important measurable quantity is the sus- 
ceptibility Xk{t,t') = J., dsRk(t,s) of the mode fc, i.e the 
total linear response to an external force of spatial modu- 
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lation k applied during the interval [t', t]. The results be- 
low imply Xk(t w + t, t w ) = k~ 2 F[k 2 xo{tw + t, t w )] where 
Xo{tw + T, tu,) -1 is a 'running' effective mass which exclu- 
sively depends on times through the local displacement, 
Xo = Xo[b], and goes to zero at large time separation r. 
The typical internal distance ~ y/xo ~ b 1 '^ grows slower 
with r as the age t w increases. 

We study fi — ► after the large-times limit. If one 
takes /x = from the start, one must take into account 
diffusion: B xx (t, 0) — + oo at large t since J dkb^ — oo. In 
addition to the aging regime defined as t, t' — > oo with 
B xx (t,t') fixed, there is then a diffusion regime where 
B xx (t,t')/B xx (t,0) is fixed Our results for S fe and 
R k are expected to hold also at fi = in the aging regime. 

Generalization of equilibrium theorems, two-time scal- 
ings. Let us describe in more detail the aging regime as 
derived below. For large times, the precise manner in 
which TTI is violated is described by 'triangle' relations 
Jl7| involving any three times: 

B k (t nl i n , t max ) f k {B k , ti n t ) , Bfc {ti n t , t max ) ) (1) 

while the violation of FDT is given by: 

R k (t,t') = X k [B k (t,t')} d t ,B k {t,t') (2) 

where (||) means that X k depends on the times only 
through B k (t,t r ). In order to complete this ansatz, we 
have to specify how f k and X k for different k are re- 
lated to one another. This is done below in an algebraic 
manner, implying no new hypotheses with respect to the 



d = case It 7 KOI. We find that 



B k (t,t') = B k [k,B (t,t')} (3) 

i.e. all fc-modes depend on times only through the depen- 
dence of one of them. Thus one can use any two-time 
function, e.g B Q (t,t') or B(t,t') to parametrise the two- 
time dependence. We also show: 

X k [B k ] = X k [B k [k,B ]] = X [B ] = X[B] (4) 

i.e. the values of all X k are the same for times such 
that B D takes the same value. We have also defined 
X[b{t,t')] = f(t,t')/d t >b(t,t'). Eq. (g) implies that if X 
is time independent in the aging regime, then it takes the 
same value for all k. Hence, for the /c-mode susceptibility 
we have Xk(t,t') = Xk[k,Xo{t,t')]. These functional de- 
pendences are testable using two-time parametric plots. 
Their explicit forms for this model are determined below. 

We find that: (i) in the FDT-regime B k < b\, 
X k [B k ] = -1/(2T) (ii) in the aging-regime b° k > 
b k > b\ there are two distinct cases. For SR corre- 
lations X k (b k ) = X and f k has the form f k (u,v) — 
3k~ l {3k{u)3h{v)). This implies B k (t,t') = j k 1 {h{t') / h{t)) 
where h(t) is increasing and independent of k. For LR 
correlations X k (b k ) is a non-constant function of b k . The 
function f k is f k (u,v) = max(u,u). 
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Formal solution for the spatial scaling. In order to 
extract the scaling properties of this model and to justify 
the above ansatz, we encode the correlation and response 
functions in the superspace order-parameter psj] . At the 
saddle point, using causality, Q xx ,(l, 2) reads 

Q xx ,{l,2) = C xx ,{ti,t 2 ) + (8 2 -6 1 ) [0 2 R xx/ (t 1 ,t 2 ) + 6 1 R xlx {t 2 ,h)] (5) 

The 6*'s are Grassmann variables and we denote 1 = 
all = dtidOiSi, 6(1 - 2) = (6» 2 - 9i)(6 2 - 
Oi)S(t 2 - h) and D<- 2 \l) = d Bl (d Sl ~ 6 x d tl ). We use 
two types of functions of the super-order parameter: 'op- 
erator' as in Q xx ,(l,2) = J d3 Q XX ,(1,3)Q XX ,(3,2) and 
'pointwise' as in Q'J(1, 2) = [Q xx ,(l, 2)] 2 . 

It is now easy to write the equation of motion for the 
order parameters, exact at the mean-field level: 



- V 2 + n + J d3V'(B xx (l, 3))) Q xx ,(l, 2) - S d (x - x')S(l - 2) - 2 [V'(B XX )Q XX ,^ (1, 2) = (6) 



where B xx ,(l, 2) = Q xx ,(l, 1) + Q xx ,(2, 2) - 2Q XX ,(1, 2). 

Using space translational-invariance one finds that all 
fourier modes Q k can be expressed in terms of the zero 
mode Qk-Q = Q through the operator relation: 

Q fc (l J 2)=[fc 2 a + 0- 1 ] _1 (l,2). (7) 

This is (in an encoded notation) a scaling relation 
for the correlation and response functions involving the 
two times and spatial separation. Substituting Eq.(^) in 
Eq.(^|) one can write a separate equation for each mode. 
Hence, we have mapped the problem of each Bf.(t,t'), 
Rk(t,t') into an effective d = problem with a compli- 
cated memory kernel f26| . Up to now we have made no 
approximations. A first important result is that one can 
now solve numerically the equation for one of the modes 
and then recover the space dependence algebraicly from 
Eq.(Q). The above property implies that an ansatz as 
in problems without space-dependence 0,|2(J applies to 
each mode B k (t,t'), R k (t,t') and justifies Eqs. 

Now, in the large time-limit any operator func- 
tion F[Q A ] of an order parameter Q A which can be 
parametrized by Ja,Xa, yields a new order parame- 
ter Q B = F[Q A ] which can also, at long times, be 
parametrized in the same form with The ex- 

plicit computation of fg, Xb in terms of /a, Xa has been 
done in Ref. pl| and when applied to Eq. (^) yields ex- 
plicit functional relations between Q k , B k and Q and, 
in particular, Eqs. (||),(||) f° r the components B k , Rk- 

One can extend this calculation to O(N) spherically 
constrained models by letting fi be a function of time 
and by imposing Q xx (l, 1) = 1 for all t\. 

Explicit calculations. We can now establish the large 
time equations in both regimes. In the FDT-regime r^ = 
— x pO T b^ (t) with xp = — 1/(2T) and one finds a single 
equation: 



5 



db F k {r) 
dr 



2T-(k 2 +M + 4x F V(b 1 ))b£(T)+4x F 4- I dr' 

« T Jo 



V\b F {r-r'))bUr') (8) 



Neglecting the time derivatives in the l.h.s. of the full dy- 
namical equations and integrating over the FDT regime 
(sec |l7],^(J) one finds the equations for the aging regime: 

= r fe (M') (k 2 + M) + V"(b(t,t')) f(t,t') +4 J ds V"(b(t >S )) f(t,s) r k (s,t') (9) 

= -b k (t, t') (k 2 + M) + ^ (V'(b x ) - V'(b(t, i')) + 4 jf (V(6(i, a)) r fc (t, s) + V"(5(t, a)) f(t, s) b k (t, «)) 

+2T-aJ ds (y\b(t,s))r k {t',s)+V''(b(t,s))f{t,s)b k (t',s))-Aj / dsV"(b(t,s)) f(t,s) b k (s,t') . (10) 

M = — 41im t ^ 00 J Q ds V"(6(i, s)) r(t, s) is the 'anomaly'. 
These equations have time-reparametrization invariance 
which prevent us from determining h(t). The 'quasi- 
static' values b , b\ follow from letting t' — > i_ in Eq.(10): 

b\ = 2T/{k 2 + M) (11) 

Similarly, Eq.(||) integrated over fc, yields either the high- 
T solution f(t,t_) = or the low-T condition 



1 = -4V"(& i ) / (k + M) . (12) 

Jk 

This hnplies 6 1 = -ATc d /(d - 2)(-4c d V" (b 1 ))^- 2 ^^-^ 
and M = (-4c d V"(6 1 )) 2 /( 4 - d ) for d < 4, where c d = 
J, (A; 2 + 1) -2 . Thus it is not necessary to know the details 
of the aging solution to determine b\, b and M, 

FDT regime. Defining 0(r) = 4(V'(6 | F (r)) - V'^ 1 )), 
the Laplace transform w.r.t. r of Eq. (pft yields: 

&£( s ) = (x F s)" 1 (/c 2 + M + s-a; F s0(s))- 1 (13) 

At small r, 6f (r) - (1 - exp(-A k r))/(xA k ) with A fc = 
k 2 + M - x F 4>{0). When b l k nt ~ l/(xA fe ) there is a 
crossover to a slower regime where one can neglect the 
term s (d/dr) and find a power law behaviour: &f (t) = 
- c ( b l) 2 T ~ & with determined by @: 

r[i-2/3]r[i-/3]- 2 =^y(6 1 ) (w) 

and F^ 1 ) = 4V"(6 1 ) 2 V"'(6 1 )- 1 %ln/ fe (fc 2 + M)- 2 . An 
explicit calculation gives ^(fr 1 ) = AlXib 1 ) in terms of 
the function X(b) defined by Eq. (|l9|). 

Aging regime in short range models. Power law mod- 
els are short range for 7 > j c — 2/(2 — d) and d < 2. 
With the ansatz X fc [6fc(i, t')] = X Eqs.(|) and © re- 
duce to a single equation for b k [t,t'). One must have 
A = -M/(4V'(6 1 )). As discussed above b fc (t,f) = 
j^ l {h{t')/h{t)). Defining u = ]nh(t) one has b k (t,t') = 
Bk{u — u') where < u < 00 and B k (0) — b\ and 
Sfc(oo) = b k . We obtain: 

= 2T - {k 2 + M + AXV'Cb x ))B k (u) + 4x4~ I du'V'(B(u - u'))B k (u r ) + 4(A - xp^VCb 1 ) - V'(B(u))) 
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Remarkably, this equation is formally similar to Eq. (g) 
though in a completely different variable. Laplace trans- 
forming w.r.t. u one gets: 

Bk(s) = bl + (Xs)- X [(fc 2 +M)- 1 - {k 2 +M - Xs^s))- 1 ] 

with 4>{u) = 4(V'(B(u)) - V'ib 1 )). 

At the beginning of the aging regime ti<l, we obtain 
B k {u) = b\- 4V"(6 1 )(fc 2 + M)~ 2 u a , and thus for t' ~ t: 

b(t, f ) - & 1 ~ \n a {h{t)/h{t')) ~ c(^)(T/t„,) a (15) 

where c(t lu ) = (ci In /),(£,„ )/d In £„,)". If /i(i) = t 5 we re- 
cover the trap- model scaling where c(t w ) is a constant 
. The exponent a is determined by: 

T[l + 2a]T[l + a]- 2 = X Yib 1 ) (16) 

For the power law model, xY{b r ) = (4 — d)7/(2(l + 7)) 
which gives a — > when 7 — » 7+ = 2/(2 — <i) and shows 
how ultrametricity in the case 7 < "f cr is approached. 

At widely separated times u — > 00, the approach to fr£ 
is described by a scaling form of k 2 u: 

where Z = J du(f>(u) is a constant, finite for 7 > 7 cr . In- 
tegrating over fc we also obtain the large time separation 
behaviour b(t,f) oc \n l - d/2 (h(t)/h(t')). 

Aging regime for long range models. These mod- 
els are solved by the ultrametric ansatz bk(t,t') = 
max(&fc(t, s), bk(s, t')), which inser ted in (|), @ leads 
pl[ to a single equation parametrized by b 

Xk(t,t')= X k[b] = (k 2 + 17 -17(b))- 1 (17) 

after simple manipulations as in pOR . We have defined 
M[6] = -4/ g 6 * V" (S')Jf [6']*' and used (§"#■ Taking a 
derivative w.r.t. 6, using c^XfclM = AT^Jc^&fc, dividing 
by X[b] and integrating over b, one gets: 

bk =bi + f bl dv . as) 

fe A (fc 2 + M-M(6')) 2 

It also implies the self-consistency condition 1 = 
4V"(6) J fe (fc 2 + M - M(b))- 2 which coincides with the 
marginality condition for the replicon. One also obtains 

X{b) = -a d V"'(fe)(-V"(&)) 2(d - 3)/(4 - d) (19) 

with ad = (4 - rf)(4c d ) 2 /( 4 - d V 2 - Tnis result, derived in 
a rather simple way, is formally identical to the result of 
the statics || and is here shown to apply directly to the 
off-equilibrium dynamics. The self-energy of the statics 
[ct](m) is thus formally identified with M — M(b). 
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One general prediction for the aging regime of mani- 
folds is the existence of a scaling form for b r (t,t') as a 
function of r and the two-times quantity b(t, t'): 

b(r,t,t') = b r=0 (t,t')H[r/b r = (t,t')^} (20) 

which holds at large scales, i.e large r and large, widely 
separated times. One has H[0) = 1. The Hartree method 
yields a remarkably simple analytical form. Using [l^ one 
finds g| d- b b k = -4V"(6)/(fc 2 + (-4V"(6) C(i )^) 2 , e.g in 
real space in d — 3 one has d^br = exp(V" (b)r / (2n)) . 

In conclusion, we described, using the Hartree approx- 
imation, the features of aging of a manifold in a random 
medium. This provides a frame of reference for the anal- 
ysis of data from experiments and simulations on realistic 
systems and will allow to determine in each case whether 
such a regime exists and up to what times. 
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